Adhesion of microcapsules 
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The adhesion of microcapsules to an attractive contact potential is studied theoretically. The 
axisymmetric shape equations are solved numerically. Beyond a universal threshold strength of the 
potential, the contact radius increases like a square root of the strength. Scaling functions for the 
corresponding amplitudes are derived as a function of the elastic parameters. 
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Introduction. ~ Microcapsules are hollow closed elas- 
tic capsules experimentally prepared as layered polyelec- 
trolyte sheets^ or through polymerization of surfactants 
coated on oil droplets immersed in aqueous solution^,. 
Their elastic constants have been measured by atomic 
force microscopj*^*^*^, deformation in shear flo'w^ or os- 
motically induced swellingSiS. For potential applications 
not only the knowledge of their elastic constants is cru- 
cial, but also, in particular, the understanding of their 
adhesive properties to other capsules or membranes, to 
substratesifliii or channel wallsi^ in microfluidic devices. 
In a first systematic experimental study of adhesion to 
glass surfaces using reflection interference contrast mi- 
croscopy, adhesion radii were measured as a function of 
capsule size and membrane thicknessi^. A systematic 
theory of such deformed shapes, however, is lacking. The 
analysis of this experiment as well as previous theoret- 
ical approaches arc adaptions of the standard textbook 
treatment which involves basically a scaling estimate of 
the deformation ignoring what are supposed to be fac- 
tors of order unitjii^. Such a simple approach predicts 
that for small deformations the adhesion radius scales as 
the square root of the strength of the potential. On the 
other hand, in the related problem of the adhesion of fluid 
vesicles dominated by curvature elasticity, a systematic 
solution based on solving numerically variational shape 
equations has shown that adhesion requires a threshold 
strength of the potentiaU^. A priori, there is no reason 
to believe that this threshold vanishes if a finite shear 
elasticity is invoked in contrast to the predictions based 
on the simple scaling picture. 




In this letter, we investigate systematically the adhe- 
sion of microcapsules in a contact potential by solving 
the corresponding shape equations. We will restrict our 
treatment to axi-symmctric shapes, see Fig. 1. We there- 
fore exclude configurations where the adhesion area buck- 
les inwards. Such configurations have been observed in 
recent simulations where, however, the range of the po- 
tential was comparable to the size of the capsulesiSiii. 
As a main qualitative result, we indeed find that a fi- 
nite universal threshold strength is required to induce 
adhesion beyond which the adhesion radius scales like a 
square root in the excess strength. The amplitude de- 
pends strongly on the elastic parameters. 

Model. - As an initial reference shape, we choose a 
sphere of radius Rq. Any deformation of this sphere costs 
an elastic energy 



(1) 



where A and fj. are the two-dimensional Lame coefficients, 
K = X + ^ is the compression modulus, and are the 
elements of the strain tensoui^^. This relation is valid 
as long as Hooke's law holds. The total energy of the 
capsule then reads 



F = + I ^ dA{2H - Cof - WA* 



(2) 



FIG. 1: Spherical microcapsule with radius Rq undergoes an 
adhesion transition to an axisymmetric shape with adhesion 
radius R* . 



The second term is the bending energy of the membrane 
where H is the local mean curvature and k the bending 
rigidity. Such a term is required to prevent sharp kinks in 
an adhesion geometry. We also introduce a spontaneous 
curvature Co = 2/Rq so that the original undeformed 
spherical capsule has no elastic energy whatsoever. Gen- 
erally a capsule can be adequately described by these 
thin shell equations as long as the shell thickness is much 
smaller than the radius. This assumption is generally ful- 
filled for capsules larger than a few microns. The third 
term is the standard adhesion energy in a contact po- 
tential of strength W where A* = ttR*"^ is the contact 
area. 

The shape equations for this model can be derived by 
first parameterizing an axisymmetric shape appropriately 
and then setting the first variation of F to zero. This pro- 
cedure can follow closely the corresponding one for fiuid 
vesicles with the additional feature that for an clastic 
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capsule tangential displacements have a physical signifi- 
cance whereas for fluid vesicles they correspond to irrel- 
evant reparameterizations of the shape and can thus be 
ignorediS. The technical details will be published else- 
where. 

The solutions of these shape equations and hence the 
"phase diagram" of adhesion depend on three dimension- 
less parameters 



KRl 
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where k and m are the scaled compression and shear 
modulus, respectively, whereas w is the scaled adhesion 
energy. In this way, all energies are referred to the scale of 
the bending energy. While typically k,m ^ 1, the scaled 
adhesion energy can easily be of order 1. Note that wc 
have not implied any volume constraint since at least for 
the polyelectrolyte shells the membrane is supposed to be 
water permeable. Introducing such a volume constraint, 
however, would not pose any fundamental complication 
but would add one more dimension to the phase diagram. 

Adhesion radius. - By solving the shape equations nu- 
merically, we find universally that the (scaled) adhesion 
radius of weakly adhered shapes behaves as 



r* = R*/Ra ~ a{k, m){w - Wc)^/^. 



(4) 



The critical strength of adhesion tw^ = 2 is independent 
of the elastic parameters and indeed the same as found 
for fluid vesicles without a volume constrainti^. For a 
weaker contact potential, the potential energetic gain of 
an adhesion disc is smaller than the cost in curvature 
energy to be paid. This balance is not modified by the 
elastic energies. 

The dependence of the amplitude a(fc,m) on the clas- 
tic parameters is best discussed in several steps. For a 
capsule without shear elasticity, m = 0, the amplitude 
a{k,0) decreases monotonically with k from the value 



ao = a(0,0) ~ 0.58 to 



a{k ^ 00,0) ~ 0.50. (5) 



For k = 0, the presence of the non-zero spontaneous cur- 
vature energy stabilizes a finite size of the contact radius. 
For fc ^ oo, the model corresponds to that of a fluid vesi- 
cle with area constraint. Both limit cases can therefore 
be checked independently using the shape equations of 
fluid vesicles without and with area constraint, respec- 
tively. Note that without shear elasticity, increasing the 
compression modulus from zero to infinity thus leads to a 
decrease in the adhesion radius of only about 15 percent 
for any given adhesion strength. 

A quantitatively more dramatic effect arises for a non- 
zero shear modulus. In Fig. 2. the amplitude a{k,m) is 
shown as a function of k for various values of m. For fixed 
k, the amplitude decreases with increasing m reaching a 
finite non-zero limit for m — > oo. Similarly for fixed m the 
amplitude also decreases with increasing k and reaches a 
non-zero value in the limit fc — > oo. The limit k oo 
can quite well be fitted by the power law 

a{k (X),m) aoo{l + m/mc)~^^^ (6) 
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FIG. 2: Amplitude a{k,m) of the adhesion radius (|1J as a 
function of the scaled compression modulus k for various val- 
ues of the scaled shear modulus m. Note that the abscissa is 
scaled logarithmically. 



with Qoc given above in cq. (jSJ and rric ~ 1.54. 

These data for general k and m almost collapse on a 
scaling plot if these clastic constants are replaced by the 
(scaled) two-dimensional Young modulus 
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and the Poisson ratio 



(7) 



(8) 



Fig. 3 reveals that the amplitude a(y, cr) is almost inde- 
pendent of the Poisson ratio as long as y 6, a condition 
fulfilled if min(A:, m) > 3. For smaller values of y, the am- 
plitude becomes the larger the smaller a. Note that all 
curves have the same limit value a(y ^ 0, c) ~ 0.58 ex- 
cept the curve for (7 = 1. The two cases k oo and 
TO = both correspond to cr = 1. In the first case the re- 
lation y = Am holds and the amplitude can be calculated 
directly from eq. In the second case, y is always zero 
and the amplitude varies from ao to floo for increasing k. 

For future reference, we extract from our data the fit 



a{y,a) 



ao 



{l + y/ycY/^ (l + y/y)" 



(9) 



cr) + 5.58(l-cr)3 and6 = 0.98 + 1.42(1- 
^. The first term is based on the in- 



withy = 1.85(1- 
cr) + 2.71(1 - cr, 
compressible case ^ discussed above with the crossover 
value yc — 6.14. The deviations for small y < 6 arc con- 
tained in the second term. This form provides a decent 
fit with a maximal deviation of 2 percent for all data 
shown in Fig. 3. 

Comparison to experimental data. - Our theory can 
now be used to reanalysc experimental data. Recently Ei- 
sner et alii^ investigated the adhesion of polyelectrolyte 
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FIG. 3: Amplitude a{y,a) for different values of u: a — ~1 
(*), cr = (o), a = 0.6 (x), cr = 0.9 (A), a = 0.99 (+), a = 1 
(•). For 2/ >, 6 the plots almost collapse on a master curve. 
For y < 6 a. stronger dependence on a arises. The continuous 
lines are fits to the data points according to eq. @. 



multilayer capsules (PMCs) on a glass surface. PMCs 
can be prepared with well-defined shell thickness, shell 
radius and surface energy. Therefore they are an ideal 
system to study the adhesive properties of microcap- 
sules experimentally. The contact area and the form of 
the adhered capsules were reconstructed by reflection in- 
terference contrast microscopy. Given the two material 
properties Young modulus and Poisson ratio the contact 
potential can be derived by fitting data sets of the con- 
tact radius vs. the capsule's wall thickness with fixed 
capsule's radius or the contact radius vs. the capsule's 
radius with fixed capsule's wall thickness, respectively. 

For thin isotropic shells of thickness h and three- 
dimensional Young modulus the two-dimensional pa- 
rameters are given by^^^ 



Y = Xs/i and 



Y,h^ 



12(1 -cr2) 



(10) 



Expressing eq. Q in these variables, we have for large 
2/ » 2/c 



R* f 
with 



{l2{l-a')y^f'^^y\w-W.y/' (11) 



Wc = 



Y,h^ 



6(l-a2)i?2 



(12) 



for the critical adhesion energy. The large y limit is ap- 
propriate since y = 12(1 — a'^)R^/h'^ ~ 10^ in these ex- 
periments. 

First, we compare our result to capsules with fixed 
shell thickness h = 25.4 nm and varying capsule's radius 
i?o and contact radius R* . In Fig. 4, the experimental 
data from Rcf. ()l3l) are shown. Using typical experimen- 
tal values for R* and i?o we can estimate w — Wc to be 
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FIG. 4: Comparison of experimental data (□, from Ref. IT3) ') 
to the fit of eq. IIH with Wc = (dashed line) and to the 
fit of eq. (IIH for finite Wc = 9 /iJ/m^ as extracted from 
experimental data using the parameters Ya = 294 MPa, h = 
25.4 nm, a — 1/3, and Rq = 10 ^im (full line). Note that the 
left part of this figure is scaled logarithmically to better show 
the influence of the critical adhesion energy Wc. 



of the order 100. Thus the critical strength of adhesion 
can entirely be neglected in cq. (|11|) . Fitting R* as a 
function of Ro with Wc = we extract for the combi- 
nation of adhesion energy. Young modulus and Poisson 
ratio the estimate W ~ ^3/(1 - ct^)^/^ • (1.4 ±0.1) • lO'^^ 
m. Estimates of the adhesion energy become possible 
using experimental values of the elastic constants. The 
Young modulus of the shell material is (294 ±30) MPa,^. 
(The value given in Ref. is severely overestimated 

and therefore not used.) The exact value of the Poisson 
ratio is unknown but it is usually between 1/3 and 1/2. 
We then get 0.39 mJ/m^ <W< 0.52 mJ/m^, choosing 
these ranges of parameters. The adhesion energies ob- 
tained this way are of the same order as the estimates in 
Ref. (0). In Fig. 4 the capsule's radius varies between 
7 and 18 fim. We can therefore estimate the critical ad- 
hesion energy Wc to be in the range from 3 ^J/m^ to 
24 fiJ/m^. In Fig. 4 both the fit with Wc = and 
the behaviour of R* vs. Rq for finite Wc are shown for 
comparison. 

Second, we compare our result with capsules with fixed 
radius i?o = 10 /um and fit the contact radius R* as a 
function of the shell thickness h (data not shown). From 
these data we obtain W ~ ^3/(1 - 0-2)1/2 . 3 ^ q g) . 
10~^2 Ya, in accordance with the previous estimate. 

Summaiy. - We have solved the shape equation nu- 
merically for elastic microcapsules with finite shear elas- 
ticity adhering to a contact potential. We have identified 
a threshold strength required for adhesion. For stronger 
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potentials the adhesion radius increases hke a square root 
with an amphtude depending on elastic constants which 
we here determined over the full range of possible param- 
eters. 
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